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Abstract 

We apply the operator product expansion to resum the leading non- 
perturbative corrections to the endpoint region of the lepton spectrum 
in inclusive semileptonic B — * X q tv decays, taking into account a fi- 
nite quark mass m q in the final state. We show that both for b — > c and 
b — > u transitions, it is consistent to describe these effects by a con- 
volution of the parton model spectrum with a fundamental light-cone 
structure function. The moments of this function are proportional to 
forward matrix elements of higher-dimension operators. The prospects 
for an extraction of the structure function from a measurement of the 
lepton spectrum are discussed. 
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1 Introduction 

The heavy quark expansion, i.e. a systematic expansion in powers of A/mg 
(we use A as a generic notation for a typical hadronic scale of the strong inter- 
actions), has by now become a widely used tool in the theoretical description 
of systems containing a heavy quark Q interacting with light degrees of free- 
dom jl], |2|, £|]. Its application to exclusive decays of heavy mesons and 
baryons has been explored already in great detail ||. Recently, the idea has 
been put forward to generalize the heavy quark expansion to obtain a QCD- 
based description of inclusive decays of heavy hadrons || [7|, || [L0| |TT|1 . 
Similarly to the case of deep-inelastic scattering, an operator product expan- 
sion is applied to the product of two local currents. The mass of the decaying 
quark provides the large momentum scale. The total decay rates can then 
be written as an expansion in inverse powers of tuq. The operators appear- 
ing at leading order have dimension three and correspond to the free-quark 
decay. The matrix elements of the dimension-four operators vanish by the 
equation of motion, and thus the leading nonperturbative corrections arise 
from dimension-five operators and are of order (A/mg) 2 . For the case of 
.B-meson decays, they can be parameterized by two low-energy parameters 
Ai and A 2 , which are related to the kinetic energy K b of the 6-quark inside 
the 5-meson, and to the mass splitting between B- and 5*-mesons. They 
are defined as |L2 



K h = ^p- = -±, m|,-m| = 4A 2 . (1) 
2m b zrrii, 

The operator product expansion has also been used to calculate differ- 
ential distributions, such as the charged-lepton energy spectrum in inclu- 
sive semileptonic decays. In this case, the relevant large momentum scale is 
Q = rrihV — q, where q denotes the momentum transfer to the lepton pair, 
and v is the velocity of the decaying hadron. After integrating over the neu- 
trino momentum, the operator product expansion is a combined expansion 
in powers of A/m^ and A/(mb — 2Eg), where Eg is the lepton energy in the 
rest frame of the decaying hadron. Over most of the available phase space 
these parameters are of similar magnitude. However, close to the endpoint, 
i.e. for (rrib — 2Eg) of order A, the second expansion parameter is of order 
unity, and a partial resummation of the expansion becomes necessary. For 
the lepton spectrum in charmless semileptonic B — > X u £ v decays, as well as 
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for the photon spectrum in the penguin-induced B — > X s 7 decays, this re- 



summation has been constructed in Refs. |L3[ |14], [15| , neglecting the mass of 
the quark in the final state. It has been shown that the leading nonperturba- 
tive effects can be related to a universal structure function, which describes 
the distribution of the light-cone residual momentum of the heavy quark in- 
side the decaying hadron. The purpose of this paper is to generalize this 
approach to the case where the mass of the quark in the final state cannot 
be neglected. For simplicity, we will work to leading order in perturbation 
theory. Radiative corrections, which have been calculated at the parton level 
in Ref. [|16|, [L7|, [LJ|, should however be included before confronting our re- 
sults with experimental data. We shall consider B — > X q t v transitions and 
treat 771 q cLS cL free parameter. As we shall see, the presence of several mass 
scales leads to technical and conceptual complications. Let us define the 
dimensionless ratio 

2 

P=-§- 2 

mi 

It will be natural to distinguish the three cases where p is of order 1, A/ra&, 
and (A/rrib) 2 , respectively. The first case is of no phenomenological impor- 
tance and will not be considered in detail. The second case is relevant for 
b — > c transitions, where p ~ 0.1. The third case applies, e.g., when one 
studies the effect of a small constituent mass of the w-quark in b — > u tran- 
sitions. 

The necessity of a resummation of the naive expansion in powers of A/rrib 
close to the endpoint is apparent from the result obtained in Refs. [|7|, ||, [| 
for the lepton spectrum in B —>■ X q i v decays. Let us define 

2E e gll^mf 
y= 7^' Tb = 192^ ' (3) 

and divide the differential decay rate into two pieces: 

The first part is the result obtained in the free-quark decay model. The 
function F(y, p) is given by 

F(y, p)=y 2 {3(1 -p)(l-R 2 )-2y(l-R 3 )}; R = ^- . (5) 
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In Fig. [I], we show F(y,p) evaluated for p = 0.08, which is an appropriate 
value for B — > X c £u transitions. The function S(y,p) in ([|) contains the 
nonperturbative corrections to a free-quark decay. The expression for this 
functions obtained by naively constructing the operator product expansion 
to next-to- leading order in A/m^ is 



S(y,p) = y 2 



Ai 



[m 6 (i - y)] 

^1 / r> 2 



(3R 2 - AR 3 



m 2 b (l-y) 
Ai 



(R 2 - 2R 3 ) - 



3A, 



m 2 b (l-y) 



(2R + 3R 2 - 5R 3 



+ 



3m 2 
A 2 



5y - 2(3 - p)R 2 + 4R- 



m 



(6 + 5y) - 12R - (9 - 5p)# 2 + 10R 3 \ Q{1 - y - p) 



+ (A/K(l-|,)]f 



(6) 



It is apparent that the operator product expansion gives an expansion in 
two parameters, A/m^ and A/(m{, — 2Ep) = A/fm^l — y)\. Over most of 
the kinematic region, the terms in S(y,p) are of order (A/m^) 2 or smaller. 
However, provided that the parameter p is of order A/m& or smaller ,[] the 
expansion becomes singular in the endpoint region: when (1 — y) is of order 
A/m,b, terms of order (A/[mb(l — y)]) n in S(y,p) become of order unity. In 
Fig. |2|, we show S(y,p) for p = 0.08, A2 = 0.12 GeV 2 , and the two cases 
Ai = —0.1 GeV 2 and A x = —0.3 GeV 2 . The low-energy parameter A 2 can 
be extracted from the known value of the B*-B mass splitting. The average 
kinetic energy, and with it the value of Ai, are not well-known, however. 
Considering the various theoretical arguments about this quantity that have 
been discussed in the literature 
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21 , we consider the two choices 



given above as reasonable "small" and "large" values for Ai- Recently, Bigi 
et al. have derived the lower bound — Ai > | A2 in a quantum mechanical 
framework 
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However, since Ai and A2 have a different dependence on 
the renormalization point, this result cannot be rigorous once short-distance 
corrections are taken into account. In fact, it is possible to construct an 



explicit counter-example to this bound using QCD sum rules \2~2\. 

1 When p is of order unity, (1 — y) cannot become small in the physical region, and 
there is no problem with (pi). 
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The effect of nonperturbative corrections is very small over most of the 
kinematic region. Close to the endpoint, however, a sharp spike of height 
~ —\i/{m h p) 2 develops. The singular behavior of the expansion becomes 
even more obvious when one goes to higher orders; at order for in- 

stance, one encounters a ^-function at y — 1 — p. Clearly, one cannot believe 
the shape of the function S(y,p) in this region. To get a more reliable de- 
scription of the spectrum, the series in A/[mfc(l — y)] has to be resummed. In 
Ref. | 13| , it has been shown that when the mass m q of the final-state quark 



is neglected, such a resummation can be performed and leads to a smooth, 
but rapidly varying, function. This so-called shape function is a genuinely 
nonperturbative object, which can be defined in terms of forward matrix el- 
ements of certain operators in the heavy quark effective theory. A complete 
resummation of the operator product expansion is, of course, a too compli- 
cated task. What can be achieved is a summation of the leading terms in 
the limit — > oo with m^il — y) kept fixed. 

In this paper we shall elaborate on this proposal and extend it to the case 
m q not being zero. Taking the limit rrtb —>■ oo with rrib(l — y) fixed, we find 
from (|) 

S(y, p) - W (3i? 2 - 4R 3 ) Q(l-y-p) + [(A/Ml - y)]) 3 } . (7) 

Our goal will be to generalize this expression to all orders in A/[m;,(i — 
y)\. In Sect. we discuss in detail the structure of the operator product 
expansion and explain the significance of three different kinematic regions, 
which we will call the free-quark decay region, the endpoint region, and 
the resonance region. Our focus here is on the endpoint region, for which we 
construct an appropriate approximation that resums the leading singularities. 
We show that it is natural to distinguish the cases where the parameter p 
in (§) is of order A/rrib or (A/int,) 2 . We emphasize that there are, in the 
context of the operator product expansion, incalculable corrections to our 
results from bound-state effects in the final state. In the two cases considered 
above, they are suppressed by a factor (A/m^) 1 / 2 or A/m^, respectively. In 
Sect. |3], we apply our formalism to inclusive semileptonic 5-meson decays 
and calculate the leading terms in the charged-lepton energy spectrum. We 
show that, in both cases, the spectrum can be written as a convolution of 
the free-quark decay rate with a fundamental light-cone structure function. 
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In Sect. f|, we illustrate our results using a simple model, which incorporates 
the known properties of the structure function and provides a description 
of the decay spectrum in terms of a single parameter. Sect. | deals with a 
discussion of possibilities to extract information about the structure function 
from experimental data on the lepton spectra in inclusive semileptonic B- 
decays. In Sect. |6], we summarize our results and give some conclusions. 

2 Resummation of the Leading Endpoint Sin- 
gularities 

By the optical theorem, any inclusive decay rate can be related to the imag- 
inary part of a transition operator T, which is defined in terms of the time- 
ordered product of two local operators. For the cases at hand, this correlator 
is of the form 

T(q 2 , vq) = -ij d 4 x e ^"^ (B(v) | T {b v (x) T, q(x), q(0) T 2 6,(0)} \B(v)) , 

(8) 

where q is the momentum transferred to the leptons, and Tj are combinations 
of Dirac matrices. The 6-quark field b v (x) is related to the conventional field 
that appears in the QCD Lagrangian by a phase redefinition: 

b v (x) = exp(—im,bV ■ x) b(x) . (9) 

This is appropriate to make explicit a trivial but strong dependence of the 
free quark field on the large mass rrib- Written in terms of the rescaled fields, 
the hadronic matrix element in (|8|) is free of large momentum scales. Hence, 
the relevant momentum transfer is 

Q = m b v-q. (10) 

Up to the difference between rrib and rriB, the variable Q 2 corresponds to the 
invariant mass of the hadronic final state. As long as Q 2 is large enough, it is 
legitimate to expand the correlator in a tower of local operators of increasing 
dimension, which are multiplied by coefficient functions that contain inverse 
powers of Q. In general, this will be an expansion in three large parameters: 
rrib, v-Q, and Q 2 . Over most of phase space, Q 2 and v ■ Q scale with the heavy 
quark mass, and the operator product expansion reduces to an expansion in 
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powers of A/mb- There is, however, a kinematic region where v ■ Q is of order 
m&, but A 2 < Q 2 < 77t|. In this case, it is consistent to work to leading 
order in A/uif,, but necessary to keep higher-order terms in Q 2 . In fact, as 
we shall see, one has to resum these terms to all orders. 

To leading order in A/mj,, the heavy quark field b v in @ can be replaced 
by the corresponding two-component spinor h v of the heavy quark effective 
theory Similarly, the physical B- meson states are replaced by the corre- 
sponding states in the mj — > oo limit. At tree- level, the leading contribution 
to the correlator is obtained by contracting the light-quark fields: 

T(q 2 ,vq) = J ^xe lQ - x {B{v)\T{h v {x)T l S q {x^)T 2 h v {{))}\B{v)) . (11) 

Here S q (x, 0) is the propagator of the g-quark in the background field of the 
light constituents of the 5-meson. The Fourier transform of S q (x, 0) has the 
form 

S q (Q) = ' 



+ lip — m q + ie 



= (9^P^m q ) {Q1 + 2 . QD _ ml _ pp+ . e) . (12) 

Note that a derivative acting on the rescaled heavy quark field corresponds to 
the residual momentum k = pb — rribV. Since the residual momentum results 
from the interactions of the heavy quark with light degrees of freedom, it is 
of order A. 

There is a lot of information that can be deduced from the structure of the 
background-field propagator. We will be interested in the region where the 
components of Q are large (of order m b ), since otherwise the operator product 
expansion breaks down. To leading order in A/m&, we can then neglect the 
covariant derivative in the numerator, as well as the term containing two 
derivatives in the denominator. This gives 

^^FW^^^' 1 ' (13) 

The term in the denominator containing the covariant derivative is of or- 
der Arrib. To see whether it is important, we have to distinguish different 
kinematic regions. When Q 2 ~ mf, the derivative term is suppressed by a 
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factor A/rrib- To leading order, the propagator corresponds to the propaga- 
tor of a free quark with momentum Q. In this region, we thus recover the 
free-quark decay model up to small nonperturbative corrections. However, 
in the endpoint region, where Q 2 — m 2 ~ Am;,, it is not possible to neglect 
the derivative term. This is why an expansion in powers of iQ ■ D/Q 2 , which 
was applied in Refs. |7|, §, |9[ to derive (|6|), becomes singular. Instead, one 
has to use (13) for the propagator in the endpoint region. For even smaller 
values Q 2 — m 2 ~ Am q , the derivative term becomes the dominant term in 
the denominator, and the operator product expansion breaks down. This 
is the resonance region, where the invariant mass Q 2 of the hadronic final 
state is close to the mass-shell of the lowest-lying resonances containing a 
g-quark. The fact that the resonance region is parametrically smaller than 
the endpoint region implies that corrections to our results, which arise from 
bound-state effects in the final state, are suppressed in the m b — > oo limit. 
This is fortunate, since these corrections matter of principle, not 

calculable using an operator product expansion. 

Let us elaborate on this last point. To see to what level one is sensitive 
to bound-state corrections in the final state, we replace the quark mass m q 
by some effective mass (m q + 5), with 5 of order A. The sensitivity to such 
effects depends upon the size of m q . For the hypothetical case m q ~ m b , 
the corrections induced by 5 ^ are of order Amj and thus of the same 
magnitude as the term containing the covariant derivative. In this case, the 
nonperturbative corrections to the free-quark decay model are incalculable. 
In the opposite limit m q ~ A, corresponding to p ~ (A/m^) 2 , we find that the 
effects corresponding to 5 ^ are always subleading. They are suppressed 
relative to the bound-state corrections in the initial state by a factor A/rrib. 
This shows that, in the analysis of the inclusive decays B — > X u lv and 



B — > X s 7 [[13|, [TJJ], the effect of a non- vanishing u- or s-quark mass, even 
as large as a constituent mass, can be neglected. A very interesting case 
in between the above is the situation where m 2 ~ Am^, corresponding to 
p ~ A/rrib- We argue that this case applies for B — > X c tv decays, where 
p ~ 0.1. Then the term in the propagator containing the derivative is of 
the same magnitude as p, and final-state corrections induced by 5 ^ are 
again subleading. However, they are only suppressed by a factor A/m q ~ 
(A/m^) 1 / 2 . We thus expect that these corrections are more important than 
in the case of decays into light (charmless) final states, and that the approach 
to the nib oo limit will be slower. 
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The above discussion shows that, for the relevant cases p ~ A/m& and 
p ~ (A/nib) 2 , the propagator (|13D can be used to resum the leading terms 
in the endpoint region. Corrections to it are suppressed in the m& — > oo 
limit. It is possible to simplify the result further. To this end, we write 
Q = (v ■ Q)(n + 5n) and choose the vectors n and 5n such that n is a null- 
vector on the forward light cone satisfying n 2 = and n ■ v = 1, and thus 
v ■ 8n — 0. In the region where the derivative term in (|13|) is important, 
i.e. for v ■ Q of order but Q 2 of order Am?,, we have n ■ 5n ~ A/m&. 
As long as we restrict ourselves to the leading term in the large-mt limit, 
we can thus neglect 5n. This leads to the final expression for the effective 
background-field propagator: 

S g (Q) ~ / dk + 5(k + - iD + ) n2 _ L ® + ™ q 2M . , (14) 
J Q + 2k + v ■ Q — m z q + xe 

where we define n-D = D + . In light-cone gauge, n-A = A + = 0, the operator 
iD + reduces to id + and corresponds to the light-cone residual momentum k + 
of the heavy quark inside the 6-meson. Note that the appearance of the null- 
vector n is closely connected to our assumption that m 2 q is of order Amj, or 
A 2 . For m 2 ~ m 2 , one would have Q 2 ~ m 2 even in the endpoint region, and 
n would be a time-like vector |15| . However, as we have argued above, we do 



not believe that this case is of phenomenological importance. From now on, 
we shall not consider it any more. 

Let us use this form of the propagator to calculate the leading contribution 
to the imaginary part of the correlator T(q 2 , q-v). We can simplify the Dirac 
structure of the hadronic matrix element by using the identity 

K r K = \ Tr (r P v ) KK- 1 - Tr ( 7m75 p v r P v ) h v 7^75 h v , (15) 

which is valid for an arbitrary matrix T. Here P v = |(1 + ^) is a projector 
onto the large components. Since the matrix element of the axial vector 
current between 5-meson states vanishes by parity invariance, we obtain for 
the leading term 



-lmT(q 2 ,vq) = -^Tr^ ($ + m q ) T 2 (1 + ^} 

x J dk + f(k+) 5{Q 2 + 2k + vQ-m 2 g ), (16) 
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where 



f(k+) = (B{v)\h v 8{k + -iD + )h v \B{v)) 



(17) 



is the leading-twist structure function that determines the probability to find 
a fr-quark with light-cone residual momentum k + inside the 5-meson [jTJj]. We 
use a mass-independent normalization of states such that (B(v) \h v h v \B(v)) = 
1. It then follows that the structure function is normalized to unity: 



dk+f(k^ 



1 . 



(18) 



It will later be convenient to introduce the Fourier transform of the structure 
function, which is given by the forward matrix element of a bilocal operator:^ 



/(f) = I dk + e-^ t f(k + ) 



{B(v)\h v {0)Pexp 



o 



h v (tn) \B{v)) . (19) 



This function obeys the normalization condition /(0) = 1. Note that the 
path-ordered exponential is absent in light-cone gauge. 

Let us recall at this point some important properties of the structure func- 



tion [13, The moments of f(k + ) are given by forward matrix elements of 
leading- twist, higher- dimension operators in the heavy quark effective theory. 
They form a set of low-energy parameters A n defined by 



A, 



dk + kl f(k+) = i n / W(0) = (B(v)\ h v (iD + r K \B(v)) , (20) 



where / (n) (0) is a short-hand notation for the n-th derivative of /(f) evaluated 
at f = 0. Using the equation of motion of the heavy quark effective theory, 
one obtains for the first three moments 



A 1 
A, 







A, 



Ai 
3 



6 



(B(v)\h v [D^g s G^] h v \B{v)) 



(21) 



2 The definition of a heavy quark structure function as the Fourier transform of the 
bilocal operator in (nw was used, in a different context, in Ref. p3). 
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where g s G^ u = i[D^, D v ] is the gluon field-strength tensor, and Ai has been 
introduced in (p. Assuming that 0.1 GeV 2 < — Ai < 0.3 GeV 2 , we find 
that 180 MeV < ^[A~ 2 < 315 MeV. This quantity is related to the charac- 



teristic width of the endpoint region |1^, |14 |. Using the equation of motion 
for the gluon field, the moment A% may be written in terms of matrix ele- 
ments of four-quark operators, which can be evaluated in the factorization 
approximation [|l|, |4|]. This leads to the rough estimate 

A 3 « -— a s f B m B « -(270 MeV) 3 , (22) 

where we have assumed fs ~ 200 MeV and a s ~ 0.4. Summarizing these 
results, we know that the light-cone structure function is centered around 
k + = 0, has a width of order 200-300 MeV determined by the average kinetic 
energy of the 6-quark inside the S-meson, and most likely (if factorization 
holds approximately) has an asymmetry towards negative values of k + . 

The support of the structure function can be deduced by observing that 
the total light-cone momentum fraction 

* = M+=^± (23) 
\jPb)+ ™b 

must be bounded between and 1. It follows that -mj < k + < rris — m&. 
Since the structure function is defined in the heavy quark effective theory, 
corresponding to the m h — > oo limit, this implies that 

- oo < k + < A , (24) 

where A denotes the asymptotic value of the mass difference between a heavy 
meson and the heavy quark that it contains, and can be identified with the 
effective mass of the light degrees of freedom interacting with the heavy quark 
[ P5fl . We expect that the support of f(k + ) for negative values of k + is slightly 
larger (because of the asymmetry) but of the same order of magnitude, i.e. 
f(k + ) should be exponentially small for k + <^ —A. 



3 Calculation of the Lepton Spectrum 

Let us now proceed to calculate the charged-lepton spectrum in semileptonic 
5-decays. The matrices Tj in (|T6| ) are of the form 7 M (1 — 75), and we obtain 
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for the leading contribution to the correlator 

- ImT^(g 2 , v ■ q) = - [ dk + f(k+) 5(Q 2 + 2k + vQ 



x 



QW + QW-g^vQ-ie^QaVp . (25) 



The next step is to express Q in terms of the lepton momentum q, and to 
contract T^ v with the leptonic tensor. A straightforward calculation leads to 
the triple-differential decay rate 



d 3 T(B ^X q iv) 
dq 2 d(v ■ q) dE e 



Gl 



V, 



gb\ 



2tt 3 



(vq- E t ) (2m b E e - q 2 ) J dk + f(k + ) 



x5 



q 2 - 2(m b + k + ) v ■ q + rri 2 + 2k + m b - m 2 , (26) 



where Ee = v -pe denotes the charged-lepton energy in the parent rest frame. 
To the order we are working, i.e. to leading order in the large-m b limit, we 
can rewrite this expression in such a way that all dependence on m b and k + 
comes through the combination 



m* b (k + ) =m b + k + , 



(27) 



which we shall identify with the effective mass of the 6-quark inside the B- 
meson. We thus observe that the leading bound-state corrections amount to 
averaging the parton-model rate for the decay of a quark with mass ml{k + ) 
over the distribution function f(k + ). The free-quark decay model is recovered 
in the limit f(k + ) — > S(k + ). 

Next we integrate over v ■ q, and then over q 2 in the kinematic region 



0<q 2 <2E t m* b (l-R*); R* 



mi{mi 



2E e ) ' 



(28) 



which follows from the requirement that < q 2 < AE^E y . We obtain 



dr Gl I V, 



qb\ 



dEp 



12tt 3 



E 2 Jdk + f(k + )& 



*2 



m 2 - 2m* b E e 



x 



{3(mf - m 2 q )(l - Rl) - Ami Ee (1 - i$} . 



(29) 



This result has several interesting properties. We first note that the heavy 
quark mass m b does no longer appear explicitly. For this reason, and in 
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particular when the focus is on the endpoint region, it would be unnatural to 
introduce the rescaled lepton energy y = 2£^/mj. Hence, we will hereafter 
present our results as function of the lepton energy E e , which is the quantity 
that is actually measured in experiments. Note, in particular, that (to the 
order we are working) the maximum value of the lepton energy is correctly 
reproduced. From the fact that &;™ ax = A = — m^ according to fl24|), it 
follows that 

2 

K) max = m B , = ^ ( 1 - ^ ) ■ (30) 

This should be compared to the kinematic endpoint Ef 1 '^ = (mb/2)(l — p) 
predicted by the free-quark decay model. The difference between the heavy 
quark mass and the physcial .B-meson mass is correctly accounted for in our 
approach. Note, however, that we are not able to account for the fact that, 
instead of the quark mass m q , there should appear in fl3HD the mass of the 
lightest meson containing the g-quark. As discussed in Sect. || this effect 
is subleading, i.e. it vanishes in the rrif, — > oo limit, whereas the difference 
between mg and rrif, remains. 

For lepton energies not too close to the endpoint, the difference between 
nib and the effective mass m£ becomes irrelevant, and (|29| ) reduces to the 
result of the free-quark decay model, as given by the first term in @. Since 
the first moment Ai in (^) vanishes, the nonperturbative corrections in this 
region are of order (A/mft) 2 . Note, in particular, that 

/ dk + f(k+) [ml(k + )} n = ml {l - ^ + . . . } , (31) 

meaning that, up to second-order corrections, it is in fact the mass of the 
6-quark that matters in the main region of phase space 0. 

It is clear from this discussion that substantial nonperturbative correc- 
tions show up only in the endpoint region, where {m^ — IEg) is of order A, and 
the difference between m^ and ml becomes important. We can separate these 
effects from the free-quark decay distribution by subtracting a term 5{k + ) 
from the structure function. In this way, we obtain the lepton spectrum as 
a sum of two terms, as shown in (|^). Introducing then the rescaled lepton 
energy y = 2Ei/m\ ) and the mass ratio p = m^/mf, and keeping only the 
leading contributions in the large-m^ limit, we obtain for the shape function 
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defined in (f|) 

A 

S(y, p)=y 2 J dk + [f(k+) - 6{k + )) (1 - 3R 2 + + less singular terms, 

(32) 

where now 

R,= ""L, ■ (33) 
m b (l -y) + k + 

This is the correct generalization of (|7]). If we would formally expand the 
structure function fllTT) as 



oo / -i \n 

f(k + ) = 5(k + ) + £ - — y-A n 6^(k + ) (34) 

n=2 n - 

and truncate the series at n = 2, we would recover the singular term shown 
in (|7|). Such a truncation is not justified, however, since every term in (|34]) 
is of the same magnitude. 

For completeness, we note that the form of the convolution (^) simplifies 
in the limit m q 0, which is relevant for B — > X u i v decays. One obtains 



dr _g 2 f \ v qb \ 



dE £ 12tt 3 
_ G F | V qb \ 



El J dk + f(k+) 6(m* - 2E e ) m* b (3m* - 4E e ) 

A 

m b E 2 e (3m b -AE e ) J dk+f(k+). (35) 



12?r 3 

2E e -m b 



This agrees with the result obtained in Refs. [13, 14, 15 



4 A Realistic Model 

At this point, it is instructive to illustrate our results with a simple, but 
realistic model. To this end, we propose the following one-parameter ansatz 
for the light-cone structure function: 

f(k + ) = —(l~x) 2 expl l --(l-x) 2 j&(l-x); z = , (36) 
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where A = ms — vn^ is treated as a free parameter. Below we shall use the 



value A = 0.57 GeV, which is predicted by QCD sum rules |26|, |27]]. Our 
model structure function is shown in Fig. |3|. It obeys all requirements that 
have been pointed out in Sect. 0. The support of f(k + ) is limited to values 
k + < A. The integral over the structure function is normalized to unity, and 
the first moment vanishes. The higher moments are proportional to powers 
of A. In particular, we obtain 

^ 2 = -y=(y-l) A 2 ~ (0.42 A) 2 , 

A 3 = -(2-y) A 3 ~-(0.33 A) 3 . (37) 

For A = 0.57 GeV, we find Ai ~ -0.17 GeV 2 and A 3 ~ -(190 MeV) 3 . These 
numbers agree well with our estimates in Sect. 0. Given these low-energy 



parameters, we can compute the 6-quark mass from the expansion QT2 

Ai + 3A 2 

»2 



m B = m b + A — -2 \- ... . (38) 



In our model, we obtain mj ~ 4.71 GeV. Similarly, we find m c ~ 1.35 GeV 
for the charm-quark mass, and hence p ~ 0.08. We shall use this set of 
parameters throughout the paper. 

In Fig. [|a, we show the lepton spectrum (l/r&) dY/dE^ obtained from (fZ9]) 
in comparison with the free-quark decay distribution. The difference between 
the two spectra is shown in Fig. ||b. It corresponds to the shape function, 
which describes the leading nonperturbative corrections to the free-quark 
decay model. For the ease of comparison with Fig. |2j, we have multiplied the 
vertical scale by mb/4 in order to comply with the definition of S(y, p) in (Q). 
In order to convert the horizontal scale to the variable y, one would have 
to multiply Ei by 2/m&. We find that the shape function is indeed sizable 
only in the endpoint region. Comparing the result of our resummation to 
the singular form of S(y, p) given in (^) and shown in Fig. we see that 
the resummation has eliminated the unrealistic spike at the endpoint. What 
remains is a smooth function, which is rapidly varying on scales of order 
Ay ~ 1. In other words, the decay probabilities have been redistributed such 
that the height of S(y, p) is strongly reduced; however, and most importantly, 
the shape function now extends beyond the parton-model endpoint (which 
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would correspond to 2.17 GeV). Needless to say, the precise shape of S(y,p) 
depends on the form of the structure function f(k+). In the following section, 
we shall discuss some strategies how to extract the structure function from 
data. 

In Fig. [5|, we show the corresponding spectra for charmless B — > X u I v 
decays, setting p = 0. The nonperturbative effects in the endpoint region be- 
come more pronounced, due to the fact that the free-quark decay distribution 
ends with a step-function in this case. 



5 Extraction of the Structure Function 

In this section, we briefly discuss how one could, in principle, extract informa- 
tion about the structure function f(k + ) from experimental data on inclusive 
semileptonic .B-decays. We should mention from the beginning that our dis- 
cussion will be incomplete in that it neglects radiative corrections as well as 
the effects of experimental uncertainties. Both could be important. It is also 
clear from the results of Refs. [|14], [T5j that a much better place to extract 
the structure function would be the photon spectrum in inclusive B — > X s 7 
decays. However, in view of the fact that this spectrum will be very hard to 
measure in the near future, whereas very detailed data for semileptonic decay 



spectra already exist ||28|| , we think it is worthwhile to consider possibilities 
to extract at least some relevant information from semileptonic decays as 
well. 

We start with a discussion of the charmless decays B — ► X u iv. They 
are particularly simple, since the nonperturbative corrections in (|35|) are 
contained in an integral over the structure function. Up to an overall nor- 
malization factor (which depends on V c b and m^), one can directly extract 
the function F(Ei) defined as 

F (Et)= m J b _ 4Et) §; (39) 

2E t -m b 

from a measurement of the lepton spectrum. The normalization can be fixed 
by observing that F(Eg) must approach unity when mj — 2Ee ^> A. Using 
the definition (|17D of the structure function, we obtain 

F(E t ) = (B{v)\h v G{m b -2E t +iD + )h v \B{v)}. (40) 
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The derivative of F{Et) with respect to the lepton energy gives the structure 
function evaluated at k + = 2Ee — m^: 

F'(E e ) = -2f(2E e -m b ). (41) 

Given a measurement of F(Eg), one can extract the moments A n of the 
structure function, which have been defined in (p0|), by integration with ap- 
propriate weight functions. One finds that [|13|] 

/ dE e (2E e - m b ) n \F(E e ) - 0(2m b - E t )] = • (42) 

J 1 J 2[n + 1) 

The fact that A\ = can be used to fix the value of m b in the step-function. 
In practice, the presence of several sources of experimental and theoretical 
uncertainties (in particular, radiative corrections and corrections of order 
A/m&, which we neglect in the above expressions) will probably limit this 
extraction method to the first few moments. 

Let us now turn to the case of B — > X c lv decays, where the form of 
the convolution integral ( p5|) is more complicated. Let us assume that the 
parameters m& and m c have been extracted from a fit of the spectrum to the 
free-quark decay model in the region far away from the endpoint. One can 
then write the observed lepton spectrum in terms of the parameter 

e = 2E e - m b . (43) 

To leading order in the large-m^ limit, the differential decay rate (p9|) can be 
rewritten in the form of a convolution of the structure function f{k + ) with 
the parton-model distribution function G(e), 

E]W e ^ F{£) = I Hk+) G{£ ~ k+) ' (M) 
where, according to (0), 

G(e) oc Q(-m b p - e) {3m 6 (l - p)(l - R 2 ) - 2{m b + e)(l - R 3 )} (45) 

up to an overall factor, and R = —m b p/e. Let us now introduce the Fourier 
transforms 

p(t) = J dee- l£t F(e), 

G{t) = f dee- t£t G(e). (46) 
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Using flT9]), we find that 

/(t) = ffil, (47) 
J 1 ; G(t) 1 ; 

i.e., by taking the ratio of the Fourier transforms of the observed lepton 
spectrum and of the parton-model distribution, one can in principle extract 
the matrix element of the bilocal operator in ([19]). The Fourier transform of 
f(t) gives the structure function. Note that the normalization of the spectra 
is irrelevant in this context, since we know the normalization of f(t) at the 
origin: /(0) = 1. 

Given an experimental determination of F(t), one can in principle com- 
pute the low-energy parameters A n using (EO), i.e. 



n (M] 
dtJ \G{t)J 



(48) 



i=0 



In practice, such an analysis will again most likely be limited to the first few 
moments. 

A less ambitious approach that could be taken is to rely on a theory- 
inspired ansatz for the structure function f(k + ), which should depend upon 
few parameters with a well-defined physical meaning. The goal would be 
to extract these parameters from a fit to experimental data. This proce- 
dure is familiar from the present analysis of inclusive decay spectra in the 
context of the phenomenological model of Altarelli et al. f29fl . We believe 



that a reasonable parameterization of f{k + ) should contain (i) the param- 
eter A, which determines the gap between the parton-model endpoint and 
the physical endpoint of the lepton spectrum, (ii) the parameter Ai, which 
is proportional to the width of the endpoint region, and (iii) an asymmetry 
parameter, which is related to the third moment A 3 of the distribution func- 
tion. An extraction of these fundamental quantities, even if it is affected by 
substantial uncertainties, would be most desirable. 



6 Conclusions 

In the endpoint region of the lepton spectrum in inclusive semileptonic decays 
of .B-mesons, a naive operator product expansion in powers of A/m& breaks 
down. The reason is that in this region of phase space, there emerges a second 
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expansion parameter, A/ (2m& — J5g), which is much larger than A/m&. In this 
case, a partial resummation of the operator product expansion is necessary 
before the theoretical results can be compared with data. This resummation 
is such that one sums all contributions of order [A/(2m& — Eg)] n , keeping 
however only the leading terms in A/m&. In many respects, it resembles the 
summation of leading-twist contributions in deep-inelastic scattering. 

For the cases of B — > X u £ v and B — > X s 7 decays, where the mass of the 
quark in the final state can be neglected, it has been shown in Refs. [|13|, [L4|, |15| 
that the leading nonperturbative contributions close to the endpoint may 
be resummed into a light-cone structure function f(k + ), which gives the 
probability to find a 6-quark with light-cone residual momentum k+ inside the 
.B-meson. This function is defined in terms of forward matrix elements in the 
heavy quark effective theory. It is thus independent of m& and has a universal 
character. In the present paper, we have shown that the same structure 
function describes the inclusive decays in the presence of a finite quark- 
mass in the final state, provided that the ratio p = {m q j m^) 2 is consistently 
treated as being of order A/mb- In this case, however, the nonperturbative 
corrections to the leading behavior are of order (A/m?,) 1 / 2 ~ A/m q . Some 
of these corrections are due to hadronization effects in the final state and 
appear to be incalculable in the context of the operator product expansion. 
Numerically, they could be quite significant. Yet, since p ~ 0.08 for case of 
B — > X c £ v decays, we argue that our results should apply and should provide 
a correct description of the leading nonperturbative effects. We should point 
out that a different position was taken by Bigi et al. [15|], where it was 



argued that b — > c transitions are close to the so-called "small velocity limit" , 
where m& — m c is assumed to be much smaller than m&, corresponding to 
p ~ 1. In this case, even the leading effects are no longer described by the 
universal light-cone distribution function. Which of the two cases is closer to 
reality remains to be seen. We note, however, that the experimental fact that 
semileptonic .B-decays are not saturated by the exclusive modes B — ► D^*H v 
indicates that the "small velocity limit", in which these two modes would 
saturate the total rate, can only be a rather crude approximation. 

In both cases, B — > X u £ v and B — > X c £v decays, the lepton spec- 
trum can be written as a convolution of the free-quark decay distribution 
with the universal structure function f(k + ). We have shown that the ef- 
fect of the momentum distribution of the heavy quark inside the meson can 
be understood in terms of an effective mass ml = rrib + k + , which deter- 
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mines the decay kinematics. Over most of phase space, the free-quark decay 
distribution is slowly varying on scales of order A, whereas the light-cone dis- 
tribution function is sharply peaked around k + ~ with an intrinsic width 
of order A. In this case, the convolution reproduces the parton model up to 
small nonperturbative corrections of order (A/m^) 2 . In the endpoint region, 
however, the free-quark decay distribution falls steeply, and the convolution 
becomes sensitive to the details of the structure function. This leads to large, 
genuinely nonperturbative effects in the lepton spectrum, which may all be 
attributed to bound-state corrections in the initial state. In particular, as 
one approaches the endpoint, the effective mass approaches the mass of the 
physical 5-meson, and one recovers the correct position of the maximum 
lepton energy. 

We have illustrated the effects of nonperturbative corrections using a sim- 
ple one-parameter model, which however includes many of the ingrediences 
of a more sophisticated description. For simplicity, and since our main fo- 
cus in this paper was to investigate the effects of bound-state corrections, 
we have not included in our discussion perturbative QCD corrections from 
the emission of real and virtual gluons. Such effects have been considered in 
Refs. [0, T!\, [18], and more recently in Refs. [|1^, 3[J. For inclusive decays into 



final states containing light quarks, their interplay with the nonperturbative 
corrections considered here is rather intricate because of the presence of large 
Sudakov double-logarithms. In particular, in B — > X u lv decays the simple 
factorization of bound-state corrections into an integral over f(k + ) [see fl35f) 1 
is replaced by a more complicated convolution of the structure function with 
a hard-scattering amplitude. For the case of B — > X c tv decays, however, 
such effects are known to be less severe. Nevertheless, further investigation of 
radiative corrections is necessary to put our formalism on a more quantitative 
basis. 

In the last part of the paper, we have discussed some possible approaches 
that could be taken to extract information about the structure function from 
semileptonic decay spectra. We are aware that such an analysis will be 
complicated, due to various theoretical and experimental limitations. It is 
likely that the most promising approach will be to make a theory-motivated 
ansatz for the structure function that contains few physical parameters, and 
to extract these parameters from a fit to data. We are confident that it 
should be possible to obtain an (approximate) determination of the first two 
non-trivial moments of f(k+). They are proportional to the average kinetic 
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energy of the 6-quark inside the 5-meson, and to the asymmetry parameter 
A3. These fundamental parameters are of sufficient interest to try such an 
analysis using the existing data on B — > X c t v decays. 

Acknowledgements: It is a pleasure to thank Daniel Wyler and Matthias 
Jamin for useful discussions. 
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Figure 1: Shape of the lepton spectrum predicted in 
free-quark decay model. We use p = (m c /mb) 2 = 0.08. 
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Figure 2: Nonperturbative corrections to the lepton spec- 
trum obtained using a naive operator product expansion. The 
solid line corresponds to Ai = —0.3 GeV 2 , the dashed one to 
Ai = -0.1 GeV 2 . 
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Figure 3: Model ansatz ( |36| ) for the structure function f{k + ), 
evaluated for A = 0.57 GeV. 
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Figure 4: (a) Charged-lepton spectrum (1/T b ) dT/dEi in 
B — > X c tv decays. The solid line is obtained from the convo- 
lution in (p9| ) using the ansatz (|36|) for the structure function. 
The dashed line shows the prediction of the free-quark decay 
model, (b) The shape function S(y, p), which is obtained from 
the difference of the two curves in (a). 
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Figure 5: Same as Fig. £|, but for B — > X u tv decays. 
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